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ABSTRACT 



Context. LS 5039 is one of the few TeV emitting X-ray binaries detected so far. The powering source of its multiwavelength emission 
can be accretion in a microquasar scenario or wind interaction in a young nonaccreting pulsar scenario. 
Aims. To present new high-resolution radio images and compare them with the expected behavior in the different scenarios. 
Methods. We analyze Very Long Baseline Array (VLB A) radio observations that provide morphological and astrometric information 
at milliarcsecond scales. 

Results. We detect a changing morphology between two images obtained five days apart. In both runs there is a core component with 
a constant flux density, and an elongated emission with a position angle (PA) that changes by 12 ± 3° between both runs. The source 
is nearly symmetric in the first run and asymmetric in the second one. The astrometric results are not conclusive. 
Conclusions. A simple and shockless microquasar scenario cannot easily explain the observed changes in morphology. An interpre- 
tation within the young nonaccreting pulsar scenario requires the inclination of the binary system to be very close to the upper limit 
imposed by the absence of X-ray eclipses. 

Key words. stars:individual: LS 5039 - X-Ray: binaries - radio continuum: stars - radiation mechanism: non-thermal 



1. Introduction 

Very high energy (VHE) gamma-ray emission in the TeV range 
has been detected in four massive X-ray binaries. PSR 1259-63 
contains a yo ung nonaccreting millisecond (ms) radio pulsar 
(lAharonian et a l. 2005a). The accreting/ejecting microquasar 
Cygnus X-1 has recently been found to be a T eV emitter 
(lAlbert et al.l2007h . LS I -i-61 303 (lAlbert et al.l2006h. suggested 
to be a fast precessing microquasar dMassi et al.ll2004]). displays 
a ch anging milliarcsecond radio morphology that Dhaw an et alj 
( 1200 6) interpreted in the context of the interaction between the 
wind of the companio n and the rel ativistic wind of a young 
nonaccreting ms pulsar (lDubusll2006h . Finally, the nature of the 
powering sou rce in LS 5039, whether ac cretion or rotation, is 
unknown (see Aharonian et al. 2005b 2006 for the TeV results). 

The system LS 5039 contains a compact object of unknown 
nature, with mass betwe en 1 .4 and 5 Mq, or biting every 3.9 days 
an ON6.5 V((f)) donor jCasares et all 1200 5). The detection of 
elongated asymmetric emission in high-resolution radio images 
obtained with the Very Long Baseline Array (VLBA) and the 
European VLBI Network (EVN) was interpreted as evidence of 
its microquasar nature, and suggested that the source was per- 
sistently producing mildly relativistic ejections with a veloc- 
ity of ~0.15c (Paredes et al. 2000, 2002). Although the X-ray 
spectra are compatible with thos e of accreting black holes dur- 
ing the so-called low/hard state dBosch-Ramon et al]|2005h . the 
radio spectra are optically thin wit h a spectral index of -0.5 
dMartfet al.|[l998t iRibo et all [19991) . Theoretical modelHng in 
the microquasar scenario has allowed researchers to reproduce 
the observed spectr al energy distributio n (SED) from radio to 
VHE gamma-rays jParedes et al.l l20()6l) . However, the lack of 



clear accretion signatures and the similarities with the SEDs of 
PSR 1259-63 and LS I -h61 303 has led other authors to model 
its multiwavelength emission using the scenario of wind interac- 
tions (Dubus 2006). One of the predictions of this kind of mod- 
elling is the periodic change in the direction and shape of the 
extended radio morphology as well as in the peak position of the 
radio core, depending on the orbital phase. 

In this paper, we report on new high-resolution radio images 
of LS 5039 in an attempt to disentangle the two possible scenar- 
ios using morphological and astrometric information. 



2. Observations 

We observed LS 5039 with the National Radio Astronomy 
Observatory (NRAO) VLBA and the Very Large Array (VLA) 
at 5 GHz frequency on 2000 June 3 and 8. We used the VLA 
in its C configuration, both as a connected interferometer and 
as a phased array. The two observing sessions, hereafter run A 
and run B, spanned from 4:30 to 12:30 UT on the corre- 
sponding dates, and were thus centered on MJD 51698.4 and 
MJD 51703.4, respectively. The orbital phases of the system 
were in the range 0.43-0.51 for run A and in the range 0.71- 
0.79 of the followin g orbital cycle for run B (using the ephemeris 
from lCasares et all2005i) . 

We performed the observations using the phase-referencing 
technique, switching between the phase reference calibrator 
J1825-1718 and LS 5039, separated 2°.41, with cycling times of 
5.5 minutes, compatible with the expected coherence time. The 
fringe finder was 3C 345. The ICRF source J 191 1-2006, located 
at 1 1°.9 from LS 5039, was observed every 22 minutes to mon- 
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itor the performance of the observations. To check the stabiHty 
of the astrometry we observed the source J1837-1532, located 
at 3°.3^ from the phase reference calibrator J1825-1718. 

We recorded the data with 2-bit sampling at 256 Mbps, at 
left-hand circular polarization. A total bandwidth of 64 MHz was 
provided by 8 sub-bands. The data were processed at the VLBA 
correlator in Socorro, using an integration time of 4 s. 

The position used hereafter for the phase reference source 
J1825-1718, obtained by means of dedicated geodetic Very 
Long Baseline Interferometry (VLBI) observations, is aj2oooo - 
18''25'"36^53237 ± ff.00021 (or +3.0 mas) and 5j2ooo.o = 
-17°18'49'.'8534 + a.'0045 (or ±4.5 mas) in the frame of ICRF- 
Ext. 1 (Craig Walker, private communication). However, this in- 
formation was not available at the time of correlation, which was 
performed for a calibrator position shifted by Aa - -1-57.0 mas 
and A6 = +21.4 mas. In the cas e of LS 5039, due to its 
proper motion (see iRibo et al.ll2002l) . the source was found to 
be Aa = +7.6 mas and A6 - -16.2 mas away from the cor- 
related phase center, located at gj20oo.o - 18*'26"' 15'!05600 and 
5j20oo.o = -14°50'54'.'2400 (from lMarti et"anil998h . 

Regarding the VLA observations, the source 3C 286 was 
used for flux density calibration, and the source Jl 825- 17 18 for 
phase calibration and to phase up the array. 

3. Data reduction 

We performed the post-correlation data reduction using the 
Astronomical Image Processing System (aips) software pack- 
age, developed and maintained by NRAO. The phased VLA po- 
sition had to be corrected by 43 cm, according to later geode- 
tic measurements. The positions of J1825-1718 and LS 5039 
were corrected using the task clcor. As recommended for phase- 
referencing experiments, we applied ionospheric and Earth 
Orientation Parameters corrections to the visibility data using 
the task clcor. A priori visibility amplitude calibration was done 
using the antenna gains and the system temperatures measured 
at each station. We then used the fringe finder to calibrate the 
instrumental phase and delay offsets. The fringe fitting (fring) 
of the residual delays and fringe rates was performed for all 
the radio sources. We found very good solutions for 3C 345, 
J 19 11 -2006 and J 1825- 17 18. Fringes for 15% and 25% of the 
baselines were missing for LS 5039 and J1837-1532, respec- 
tively. Typical data inspection and flagging were performed. The 
obtained self-calibrated images of J 19 11-2006, which show a 
one-sided jet, are very similar to those present in the VLBA cali- 
brator source catalog, assessing the reliability of the observations 
and the initial data reduction. 

For each run, we performed standard phase and amplitude 
calibration and phase self-calibration steps on LS 5039. As a 
compromise between angular resolution and sensitivity, we used 
a weighting scheme with robust to produce the images with 
the task imagr. The obtained image had a synthesized beam of 
~6.4x2.2 mas in position angle (PA) of ~ 0°. Extended emis- 
sion was visible for both runs. To enhance the small-scale mor- 
phology, we produced final images using a convolving beam of 
3.4x1.2 mas in PA of 0° (similar to what would be obtained 
by using robust -2 within aips or a uniform weighting scheme 
within DiF MAp). This con volving beam is the one we used pre- 
viously in iParedes et alj (l2000i) . allowing for a straightforward 
comparison of the images. 

To measure the position of LS 5039 in both runs we used 
the phase-referencing technique, and transferred the phases of 
the calibrator J1825-1718, fitted with the task fring, to our tar- 
get source. The peak flux density of this extragalactic phase- 



reference calibrator is 118 mJy beam ' at 5 GHz, enough to 
provide useful phase information. However, it suffers scatter- 
ing and becomes resolved with the VLBA, providing relatively 
noisy phase information on the longest baselines. We inspected 
and screened the phase delays and delay rates of J 1825- 17 18. 
Tropospheric corrections were unsuccessful because the obser- 
vation strategy was not adequate to account for them. Phase- 
referenced images of the target source with a synthesized beam 
of ~4.9xl.5 mas in PA ~ 0° were obtained, and the peak position 
was measured with the aips task jmfit. 

4. Results 

4.1. VLA 

The VLA data of LS 5039 were compatible with a point-like 
source for the obtained synthesized beam of 5.5 x 3.7" in PA of 
4°. We measured the flux density of the source every 30 min- 
utes and obtained a mean of 29.4 mJy with a standard deviation 
of cr= 1.1 mJy during the 8 hours of run A, and 28.4 mJy with 
cr=0.7 mJy for run B. Therefore, LS 5039 displayed flux den- 
sity variability below 4% at the Icr level within each 8-hour run. 
Similar results were obtained within 4-hour halves. This guar- 
antees that we can produce reliable VLBA images for each of 
the entire observation periods, as well as for 4-hour intervals, 
provided the morphology does not change within such intervals. 

4.2. VLBA 

We show the final VLBA+phased VLA self-calibrated images 
in Fig. [1] The image obtained for run A displays a central core 
and bipolar and nearly symmetric extended emission with PA 
^ 116 + 2°, with enhanced emission toward the southeast. The 
total flux density recovered by the VLBA, obtained with the task 
TVSTAT within aips, is ~25 mJy, representing ~85% of the VLA 
value. The peak flux density of the core is 10.5 mJy beam 
The image is similar to the one obtained with the same ar- 
ray in 1999 May 8, corresponding to orbital phases 0.12-0.15, 
which showed a slightly mor e asymmetri c extended emission in 
PA~125°, as can be seen in Paredes et al.l (|2000). In contrast, the 
image obtained for run B displays a core and bipolar but clearly 
asymmetric structure, with PA ^ 128 ± 2° and enhanced emis- 
sion toward the northwest. The total flux density recovered is 
~24 mJy, or ~85% of the VLA value. The peak flux density of 
the core is 10.5 mJy beam"', as in run A. 

To characterize the extended emission we used uvfit and jm- 
fit within AIPS, as well as model fitting tools within difmap to 
check the reliability of the obtained results. The preferred model 
for run A data consisted of three Gaussian components to ac- 
count for the core (Corel), the southeast (SEl) and the north- 
west (NWl) components. A similar model could fit the data for 
run B, although the southeast (SE2) component is marginally fit- 
ted in this case. The fitted parameters are quoted in Table [T] As 
can be seen, the peak flux density of each component is ~50% of 
its total flux density, although this is partially due to the convolv- 
ing beam size being smaller than the synthesized beam size. In 
particular, we stress that the core components are not point-like. 

Finally, we splitted the 8 hours of data in each run into 4- 
hour data sets. For each of them, we performed a similar cal- 
ibration process as the one described above, and we modelled 
the final visibilities of the self-calibrated data. We measured no 
significant morphological differences between the two halves in 
any of the two runs. In particular, the distance between the fitted 
components Corel and SEl is stable in 4 hours within the errors 
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Fig. 1. VLBA+phased VLA self-calibrated images of LS 5039 at 5 GHz obtained on 2000 June 3 (left) and 8 (right). North is up and 
east is to the left. Axes units are mas, and the (0,0) position corresponds to the source peak in each image. The convolving beam, 
plotted in the lower left corner, has a size of 3.4x1.2 mas in PA of 0°. The first contour corresponds to 5 times the rms noise of 
the image (0.08 and 0. 11 mJy beam"' for run A and B, respectively), while consecutive ones scale with 2'^^. The dates and orbital 
phases are quoted in the images. There is extended radio emission that appears nearly symmetric for run A and clearly asymmetric 
for run B, with a small change of ~10° in its position angle. We have labeled the components fitted to the data, the parameters of 
which are listed in Table [T] 

Table 1. Parameters of the Gaussian components fitted to the data for each 8-hour run. Columns 3 and 4 list the peak and integrated 
flux densities of each component. Columns 5 to 8 list the polar and Cartesian coordinates of the components with respect to the 
peak position. Columns 9 to 11 list the sizes of the Gaussian major and minor axes and the PA of the major axis. The parameters PA 
and PAAxis are positive from north to east. 



Run 


Comp. 


Peak 55 ghz 


S 5 GHz 


r 


PA 


Aa 


AS 


Maj. Axis 


Min. Axis 


PAAxi. 






[mJy beam"'] 


[mJy] 


[mas] 


n 


[mas] 


[mas] 


[mas] 


[mas] 


["] 


A 


Corel 


10.54 + 0.08 


20.0 ± 0.2 










3.69 ± 0.03 


2.09 ± 0.02 


103 ± 1 




SEl 


1.11 +0.08 


2.6 ±0.2 


3.67 ± 0.08 


115.9 ± 1.7 


3.30 ±0.07 


-1.60 ±0.12 


4.1 ±0.3 


2.36 ±0.16 


17 ±5 




NWl 


0.88 ± 0.08 


1.5 ±0.2 


3.29 ± 0.09 


-63 ±2 


-2.92 ± 0.08 


1.52 ±0.14 


3.6 ±0.3 


1.99 ±0.18 


3±6 


B 


Core2 


10.45 ±0.11 


17.6 ±0.3 










3.71 ± 0.04 


1.8 ±0.2 


180 ± 1 




SE2 


0.75 ±0.11 


1.8 ±0.4 


2.8 ±0.2 


129 ±5 


2.17 ± 0.13 


-1.8 ±0.3 


4.6 ± 0.7 


2.1 ±0.3 


1 ± 7 




NW2 


2.22 ±0.11 


3.9 ±0.3 


2.94 ± 0.06 


-52.2 ± 1.4 


-2.32 ± 0.04 


1.80 ±0.09 


4.3 ± 0.2 


1.68 ±0.08 


174 ± 2 



(cr^ = 0.31 mas, cr^ = 0.62 mas). This is also the case between 
Core2 and NW2 components {cr^ - 0.55 mas, cr^ = 0.25 mas). 

4.3. Astrometry 

The phase-referenced image of LS 5039 obtained for run A, 
not shown here, displays a compact source with peak flux den- 
sity of 6.7 mJy beam"', and a faint extended emission toward 
the southwest (~2 mJy beam"'; the rms is 0.3 mJy beam"'). 
The measured position of LS 5039 for run A is ffjioooo - 
18''26'"15^05653 + 0!00001 (or ±0.15 mas), and dnooo.o = 
-14°50'54'.'2564 + 0'.'0015 (or +1.5 mas) (see discussion below). 
We recall that the reference coordinates are those ofJ1825-1718 
listed in Sect.|2] which have systematic uncertainties of +3.0 mas 
and +4.5 mas, respectively. For checking purposes we splitted 
the data in blocks of 4- and 1-hour lengths, made images for 
each block, and measured the corresponding positions. The peak 



position of LS 5039 appears to move Aa — +0. 1+0.1 mas, 
A6 = -1-2.8 + 0.2 mas between the two 4-hour blocks, and slightly 
more in the 1-hour blocks. The observed direction has a PA of 
1.7 + 1.7°, which is the same as the line joining the positions of 
the source and the phase-reference calibrator J 1825- 17 18. On 
the other hand, the expected error for differential astrometry is 
given by the separation d in degrees from the phase-reference 
source and the offset A of its correlat ed position, according to 
A X (t//180) X 7T (see IWalkerlll999l) . Plugging our off^set of 
60.9 mas and a distance of 2°.47 we obtain an error of 2.6 mas, 
very similar to the observed displacement of LS 5039. A similar 
procedure using 4-hour blocks for the astrometric check source 
J1837-1532 (d = 3^38) reveals motion of the peak of this ex- 
tended source by 3.2 ± 0.4 mas at PA=46 + 8°, as expected. 
Therefore, these secular motions appear to be purely instrumen- 
tal. The errors assigned to the coordinates of LS 5039 quoted 
above are half of the total secular motion measured in 1-hour 
blocks. 
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The phase-referenced image of run B, not shown here, re- 
veals a double source in the eastwest direction. The two compo- 
nents have similar peak flux densities of 3.6 + 0.1 mJy beam"', 
much lower than the one of the self-calibrated core, and are sep- 
arated 2.5 + 0.1 mas in PA of 87 + 5°. Moreover, 4- and 1-hour 
blocks reveal a fading of the western component and a brighten- 
ing of the eastern one along the run (as well as a similar secular 
motion as in run A). This symmetric double structure is in con- 
trast to the asymmetric structure seen in the self-calibrated im- 
age shown in Fig.[T}right. Tropospheric errors, which affect the 
phase-referenced image and cannot be accounted for, can eas- 
ily split Core2 into the observed double source. Therefore, the 
precise position of the peak of LS 5039 cannot be measured in 
run B. 

5. Discussion 

The observations of LS 5039 reported here, obtained with the 
VLBA on two runs separated by 5 days, show a changing mor- 
phology at mas scales. In both runs there is a core component 
with a constant flux density within errors, and elongated emis- 
sion with a PA that changes by 12 + 3° between both runs. The 
source is nearly symmetric in run A and asymmetric in run B 
(see Fig.[TJ- 

In the microquasar scenario, and assuming ballistic motions 
of adiabatically expanding plasma clouds without shocks (see 
iMirabel & Rodrfguezlll999l) . the morphology of run A can be 
interpreted as a double-sided jet emanating from a central core 
with the southeast component as the approaching one. The rela- 
tive distances of the components to the core would imply a bulk 
motion with velocity (3 > 0.06 + 0.02, whereas the flux asymme- 
try implies /3 > 0.08 + 0.02 for discrete ejections. The angle 
between the jet and the line of sight is restricted to <87°. These 
results are similar to those previously found with VLBA and 
EVN images dParedes et alj|2000l |2002|) . We can compute, for 
different angles, the expected displacement of the SEl compo- 
nent with respect to Corel in 4 hours using the measured value 
for the flux asymmetry of yScos - 0.08, and compare it to the 
upper limits for this displacement quoted at the end of Sect. 14.21 
This provides a constraint of 6 < 48° at the 3cr level. In contrast, 
in run B the northwest component appears to be the approach- 
ing one emanating from the core. The distances from Core2 to 
the components NW2 and SE2 are very similar and do not imply 
any significant relativistic motion. The flux asymmetry implies 
yS > 0.11 + 0.03 and < 84° for discrete ejections. As for run A, 
our measured upper limits of any displacement in 4 hours imply 
< 45° at the 3cr level. In conclusion, the lack of proper mo- 
tions implies that, for the measured flux asymmetries, the SEl 
jet should be pointing at < 48°, and the NW2 jet at < 45°. 
In this context, jet precession is needed to explain this behav- 
ior. If the precession axis is close to the plane of the sky, as in 
SS 433 (Blundell & Bowler 2004), the precession cone should 
have a semi-opening angle > 45° to fulfill the limits quoted 
above. If the precession axis is close to the line of sight, a small 
precession of a few degrees could explain the images of runs A 
and B. However, in both cases the PA of the jet should vary con- 
siderably, in contrast to the small range cov ered by all observed 
value s at mas scale, between 1 15 and 140° jParedes et al.ll2000l 
l2002h . 

Alternatively, the morphology detected in run B could be the 
result of a discrete ejection with core suppression where Core2 
is the approaching component and NW2 the receding one, while 
there is no radio emission at the origin of the ejection. In this 
case, the flux asymmetry implies /3 > 0.21 +0.01 and < 78 + 1°. 



The measured upper limits of the relative motion of these com- 
ponents in 4 hours provide a strong constraint of < 14° at 
the 3cr level. Imposing cos = 0.21, the origin of the discrete 
ejection (i.e., the position of the binary system or the position 
of Corel in run A) should be placed at Aa = -1.4 mas and 
A6 - +1.1 mas from Core2 in Fig[T] Although this would be 
self-consistent if the precise position of Core2 corresponds to 
the east component of the phase-referenced image of run B, it 
could be a chance coincidence produced by tropospheric errors. 
One should also invoke jet precession on short timescales to ex- 
plain the change in PA of the extended emission. However, large 
X-ray and radio flux density variations are observed in micro- 
quasars displaying discrete ejections (see Fender 2006, and ref- 
erences therein), while the peak and total ra dio flux densities 
of LS 5039 are s trikingly constant (see also. iRibo et al.l [r999l : 
IClark et 3111200 lb . and there is no evidence of an X-ray flare in 
11. 5 years of RXTE/ ASM data. 

In the young nonaccreting pulsar scenario, the different mor- 
phologies we have detected at different orbital phases could be 
due to the change of the relative positions between the pulsar and 
the companion star along the orbit (see Dubus 2006 for details). 
However, observations at different orbital phases have always re- 
vealed a very similar PA for the extended emission, which can 
only be observed if the binary system is seen nearly edge on. On 
the other hand, the absence of X-ray eclipses places an upper 
limit of / $75° in this scenario (Dubus 2006). Therefore, these 
two restrictions imply an inclination angle that should be close 
to the upper limit of 75°. We note that this is much higher than 
the inclination angle if the binary system is pseudo-synchronized 
dCasares et al.ll2005h . 

In conclusion, a simple and shockless microquasar scenario 
cannot easily explain the observed changes in morphology. On 
the other hand, an interpretation within the young nonaccreting 
pulsar scenario requires the inclination of the binary system to 
be very close to the upper limit imposed by the absence of X-ray 
eclipses. Precise phase-referenced VLBI observations covering 
a whole orbital cycle are necessary to trace possible periodic 
displacements of the peak position, expected in this last scenario, 
and to obtain morphological information along the orbit. These 
will ultimately reveal the nature of the powering source in this 
gamma-ray binary. 
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